Astronomy Sz Astrophysics manuscript no. bbhyp'acc 


February 2, 2008 


(DOI: will be inserted by hand later) 





Spin-up of the hyperon-softened accreting neutron stars 

J.L. Zdunik, P. Haensel, and M. Bejger 

N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warszawa, Poland 
e-mail: haensel@camk.edu.pl, jlz@camk.edu.pl, bejger@camk.edu.pl 



in 
o 
o 

(N 
£ 

O 

(N 

(N 
> 
O 

(N 

o : 

in , 
o ■ 

^ : 

Oh' 
I 1 

o 1 

a: 



Received xxx Accepted xxx 

Abstract. We study the spin-up of the accreting neutron stars with a realistic hyperon-softened equation of 
state. Using precise 2-D calculations we study the evolutionary tracks of accreting neutron stars in the angular- 
momentum - frequency plane. In contrast to the case of spinning-down solitary radio-pulsars, where a strong 
back-bending behavior has been observed, we do not see back-bending phenomenon in the accretion-powered 
spinning-up case. We conclude that in the case of accretion-driven spin-up the back-bending is strongly suppressed 
by the mass-increase effect accompanying the angular-momentum increase. 
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1. Introduction 

Some theories of dense matter predict a phase transi- 
tion to a new "exotic" state of matter at densities higher 
than the nuclear saturation density. A phase transition 
softens the equation of state (EOS). It has in turn been 
suggested that it leads to characteristic astrophysical sig- 
natures, for example, very large values of the braking 
index of spinning-down radio-pulsars, as well as period- 
clustering of spinning-up neutron stars in low-mass X-ray 
binaries (Glendenning et al. 1997, Glendenning & Weber 
2001a, Glendenning & Weber 2001b, Poghosyan et al. 
2001, Glendenning and Weber 2002). Detection of any of 
these specific features might serve as a clue to a deeper 
understanding of the EOS of matter at supranuclear den- 
sities. 

This work is a consecutive step in our studies of the re- 
lation between the realistic EOS of dense matter, dynam- 
ics of neutron stars, their evolutionary tracks, and obser- 
vational signatures. As in our previous article (Zdunik et 
al. 2004, hereafter referred to as Paper I) we would like to 
focus on the phenomenon of back-bending, which is a very 
characteristic consequence of the softening of the EOS on 
the spin evolution of a neutron star. 

In Paper I we showed that the back-bending in 
spinning-down solitary pulsar is not necessarily evidence 
of an "exotic" phase of dense matter. Specifically, it was 
demonstrated that the occurrence of back-bending can 
also be caused by hyperons, particles whose presence in 
dense matter was predicted over forty years ago (Cameron 
1959). We also demonstrated the crucial importance of 
precision in numerical calculations and of checking the 
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Fig. 1. Solid lines NH: pressure p vs. mass density p for 
the EoSl NAE and EoS2 NAE of Balberg et al. (1999), 
calculated by Balberg and Gal (1997) ( N1H1 and N2H1 
EOS in the notation of Paper I). Dotted lines were ob- 
tained by suppressing the hyperons, i.e., by assuming 
purely nucleon dense matter. The visible softening around 
p ~ 7 x 10 14 for the NH models results from the appear- 
ance of hyperons in dense matter. 



stability of quasi-stationary rotating configurations in the 
back-bending region. 
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In the present paper we describe numerical simulations 
of the mass and angular momentum transfer onto a neu- 
tron star with the EOSs involving hyperons. In particu- 
lar, we consider the equation EoS2 NAE (N2H1 in nota- 
tion of Paper I) of Balberg et al. (1999), which according 
to Paper I, revealed the strongest back-bending effect in 
spinning-down neutron stars. We present results obtained 
for two EOSs softened by the appearance of hyperons (out 
of four EOSs calculated by Ba lberg et al. 1999| | which cor- 
respond to two different models of purely nucleon matter. 
One of them, Nl, has a standard incompressibility of nu- 
clear matter K = 240 MeV. The second model, N2, corre- 
sponds to a rather stiff nuclear matter with K = 320 MeV. 
Consequently, N2 gives an equation of state that is non- 
causal at very high density. However, it should be stressed 
that all stellar configurations presented in our paper have 
central density well below this non-causality threshold. 

Our model is described in Sect. |21 Results of numerical 
calculations are collected in Sect. [31 while Sect. 01 contains 
discussion and concluding remarks. 

2. Astrophysical model 

We assume that the evolution of an accreting neutron star 
can be represented as a sequence of stationary rotating 
configurations of increasing baryon mass. In what follows 
we assume that matter is accreted from an accretion disc 
in which matter elements move on Keplerian orbits, mov- 
ing slowly inward due to the viscous angular momentum 
outward transport. The magnetic field is assumed to be 
sufficiently weak so as not to affect the disc dynamics. 
However, it may still play an important role in the an- 
gular momentum transport. The inner part of the disc is 
assumed to be geometrically thin, so that the approxima- 
tion of accretion of matter in the equatorial plane is good. 

Two radii of co-axial circles in the equatorial plane play 
a crucial role in our calculation. The first one, denoted 
i? oq , is the radial coordinate of the neutron star surface 
at the equator. The second important radial coordinate is 
the radius of the marginally stable orbit of a particle in 
the space-time around and within a rotating neutron star, 
r ms . Orbital motion of particles on the orbits with r > r ms 
is stable while that of particles on circular orbits with 
T < r ms is unstable. Two accretion regimes are considered, 
as explained below. 

The first scenario corresponds to r ms > R eq . We then 
calculate the spin-up of an accreting neutron star accord- 
ing to the prescription given by Zdunik et al. (2002). 
Matter is accreted from the marginally stable orbit (MSO) 
with r = r ms , which determines the inner edge of the ac- 
cretion disc. A particle falls from the MSO onto the neu- 
tron star surface after losing an infinitesimal amount of 
its angular momentum due to viscous processes. We de- 
note the specific angular momentum per unit baryon mass 
at the MSO by l ms . A more detailed discussion concern- 
ing determination of MSO and l ms in the framework of 
General Relativity can be found in Zdunik et al. (2002). 
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Fig. 2. Gravitational mass versus central baryon num- 
ber density along several evolutionary tracks of accreting 
stars with hyperons (model N1H1) and without hyperons 
(model Nl). Solid lines (model N1H1) and dotted lines 
(model Nl) correspond to the typical case - all angular 
momentum of matter falling from the marginally stable 
orbit is transferred to the star (xi = 1). Dash-dotted lines 
represent three cases for an Mq = 1.46 M Q N1H1 star with 
xi ~ 0.9,0.7,0.5, fractions indicated as boldface labels, 
from the left to the right. The long-dashed, nearly vertical 
line at maximum mass is defined by the onset of instabil- 
ity with respect to axi-symmetric perturbations. Dashed 
lines correspond to the condition R cq — r mB . Termination 
points are marked by a filled circle for Keplerian configura- 
tions and by an open circle for the instability with respect 
to axi-symmetric perturbation. 

In the second case r ms < R cq , and MSO is located in- 
side the star. In this case we assume that the inner edge of 
the accretion disc is connected with the stellar surface via 
a thin transition layer, where viscous processes transport 
the angular momentum to the neutron star. 

In both cases we define the innermost stable circu- 
lar orbit (ISCO), from which the angular momentum 
and baryon mass is transferred to a neutron star. The 
ISCO determines the inner edge of the accretion disc. We 
use the approximation nsco = r ms for r ms > R cq and 
nsco = -Rcq for r ms < R cq . 

The value of specific angular momentum per unit 
baryon mass of a particle orbiting the neutron star at the 
ISCO, /is, can be calculated by solving exact equations of 
the orbital motion of a particle in the space-time produced 
by a rotating neutron star, given in Appendix A of Zdunik 
et al. (2002). 

Accretion of an infinitesimal amount of baryon mass 
(1Mb onto a rotating neutron star is assumed to lead to 
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Fig. 3. Angular momentum as a function of spin fre- 
quency for accreting, spinning- up NS, with xj = 1 for 
EOS with (solid lines, N1H1 model) and without hyper- 
ons (dotted lines, Nl model). The label gives the baryon 
mass for the starting non-rotating configuration. The J(f) 
track of an isolated NS with baryon mass of 1.8 Mq (thin 
dashed line) is also shown. To make comparison easier, 
the initial rotating state of spinning down isolated NS has 
been chosen to lie on the spin-up track of accreting NS 
with initial baryon mass of 1.46 Mq. Termination points 
are marked by a filled circle for Keplerian configurations 
and by an open circle for the instability with respect to 
axi-symmetric perturbation. 

a new quasi stationary rigidly rotating configuration of 
mass Mb + dMe and angular momentum J + d J, with 

dJ = asjlia dAf B , (1) 

where xi denotes the fraction of the angular momentum of 
the matter element transferred to the star. The remaining 
fraction 1 — xi might be lost via radiation or other dissi- 
pative processes. In the case of r ms > i? cq , these processes 
are acting on the MSO and during the fall of a matter 
element through the dilute but magnetized plasma and 
the photon radiation emitted from stellar surface, which 
fills the gap between the MSO and stellar surface (Agol 
& Krolik 2000, Afshordi & Paczynski 2003). In the case 
of nsco = Rcq, the dissipation and the angular momen- 
tum loss take place within the boundary layer, through 
which the accretion disk joins the stellar surface (see, e.g., 
Popham & Sunyaev 2001, and references therein). 

3. Calculations of evolutionary tracks 

We assume that the evolution of an accretiting neutron 
star can be represented as a sequence of stationary ro- 
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Fig. 4. Gravitational mass versus the (circumferential) 
equatorial radius, along several evolutionary tracks of ac- 
creting strange stars for N2H1 model. The dotted line rep- 
resents the model for which only 75% of angular momen- 
tum is transferred to the star. The long-dashed line at 
maximum mass is defined by the onset of instability with 
respect to axi-symmetric perturbations. The dash-dot line 
corresponds to the Keplerian (mass shedding) limit. The 
dashed line at R cq ~ 13 km corresponds to the condition 



fating 2-D configurations. The numerical method used for 
the calculation of configurations belonging to this family 
was described in Sect. 3 of Paper I. Here, we would like 
to stress the very high precision of our calculations; eval- 
uation of the GRV2 and GRV3 virial error indicators (see 
Nozawa et al. 1998 and references therein) showed relative 
errors within ~ 10~ 5 . 

To study the effect of hyperons on the evolution of ac- 
creting neutron star, we compare the results for Nl and 
N1H1 EOSs, in the notation of Paper I, models of Balberg 
& Gal (1997). In Fig.FJlwe present the evolutionary tracks 
in n c —M plane, where n c is the central baryon density and 
M is gravitational mass of NS. Using this figure we can 
answer the question how accretion changes physical con- 
ditions in the center of a star. As we can see for xi = 1, 
the curves for N (nucleons only) stars are almost vertical. 
It means that accretion results in the increase of / (and 
J), but the density in the center of the star remains essen- 
tially unchanged. In the case of NH EOS (with hyperons), 
at the first stage of accretion the density in the center in- 
creases, as well as the size of the hyperon core. However, 
the main effect of the mass and angular momentum trans- 
fer is spin- up and change of shape (deformation) of NS, 
accompanied by a decrease of n c . As in the case of strange 
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Fig. 5. Spin frequency as a function of angular momen- 
tum of the solitary, spinning-down pulsars (dashed lines) 
and accreting, spinning- up NS with xi — 1 (solid lines). In 
the case of dashed lines, the numbers denote their baryon 
mass Mb in the units of Mq, which stays constant along 
the evolutionary track. For accreting neutron-star spin- 
up tracks, the label gives the baryon mass for the starting 
non-rotating configuration. The arrows indicate the evolu- 
tion of neutron-star configurations with time. Termination 
points are marked by a filled circle for Keplerian configura- 
tions and by an open circle for the instability with respect 
to axi-symmetric perturbation. For the 2.1M Q curve, no- 
tice a tiny thick segment, corresponding to the relativistic 
spin up by the angular momentum loss, just before the 
axi-symmetric instability limit is reached. Note the change 
of sign of d//dJ for the dashed (isolated pulsars) curves, 
which signals the end point of the back-bending stage, 
denoted by thick dashed line segments. There is no back- 
bending for accreting stars of similar baryon masses, for 
which d//dJ stays positive along all evolutionary track. 
For comparison we plotted one curve with xi = 0.75 (dot- 
ted line). 



stars (Zdunik et al. 2002), the effect strongly depends on 
the fraction, xi, of the angular momentum transferred to 
the stars from the ISCO. For xi — 0.5, central density and 
the size of hyperon core increase during accretion, lead- 
ing eventually to a collapse of the NS into a black hole; 
this is to be contrasted with the xi > 0.7 cases, where the 
termination point is the mass-shedding (Keplerian) limit. 
Considered examples show that, in order to build a rather 
large hyperonic core, a significant amount of angular mo- 
mentum has to be lost in the course of accretion. 

Another difference between the evolution of accreting 
NS with and without hyperon core is visualized in Fig. 



For the Nf and NIHI EOSs, the evolution of spin fre- 
quency, as a function of angular momentum, depends on 
the mass of initial configuration in opposite ways. For 
-Mb, in = I.46M , the evolution for these two EOSs does 
not differ significantly, and the curves almost coincide. 
The Nl stars with higher initial mass have larger J for 
the same /. If a star contains a hyperon core, as for the 
N1H1 model, higher mass results in a lower J for a given 
spin frequency. This effect can be explained by the specific 
dependence of the moment of inertia / on the mass of the 
star. For slow rotation, the moment of inertia of the star 

is defined as a limit / = lim dJ/df. Let us recall that, 

/-o 

except for M close to maximum allowable mass for non- 
rotating NSs, M^J 1 *, the value of / increases monotonically 
with M. However, close to M^*, a rapid decrease of R 
with increasing M dominates and makes / decrease with 
increasing M. For the Nl model the initial configuration 
with Mb ~ 1.5 — 2Mq is sufficiently far from the maxi- 
mum mass and the increase with M dominates. For the 
N1H1 model we are close to M max and the decrease of R 
with a relatively small change (increase) in M dominates. 

In Fig. |21 we also plot the sequence of spinning down 
configurations with hyperon core for isolated neutron star 
(fixed baryon mass Mb). To make comparison easier, the 
baryon mass for this sequence (1.8 Mq) is equal to the 
final mass of the accreting model with the initial mass 
1.46 Mq. The difference between the J(f) dependence is 
evident; for isolated NS during slowing down the pulsar 
loses about half of its initial angular momentum without 
changing its rotation period much (in this case ~ 800 Hz). 
This effect is almost completely suppressed by accretion, 
as the derivative d J/ df does not change significantly along 
the curve M B , in = 1.46M @ 

The same effect can be seen for the N2H1 EOS, for 
which we obtained the most pronounced back-bending 
behavior for the isolated, spinning down neutron stars 
(Paper I). Results are presented in figures 0] and [5] In 
Fig. 0] we also plot the curve defined by the condition 
Req = ^ms, which is almost vertical and corresponds to 
the equatorial radius 13 km. Marginally stable orbit ex- 
ists to the left of this curve (i.e. equatorial radius smaller 
than ~ 13 km). The evolution of accreting NS, when 75% 
of the angular momentum is transferred to the star, can 
lead to the maximum accreted mass larger by more than 
0.1M Q than in the case of xi = 1 (typically the mass in- 
crease is ~ 0.4Mq), before the star reaches the Keplerian 
limit (see Fig. In Fig. |S|we use angular momentum 
as an argument (abscissa), because the change of angular 
momentum is directly connected with a time evolution, by 
Eq. in the case of accretion and by the equation de- 
scribing the energy and angular momentum loss due to the 
physical processes (e.g. magnetic dipole radiation) for iso- 
lated spinning down NS. The value of hs does not change 
significantly during accretion period, so that the change 
of J is roughly proportional to the accreted mass (i.e., to 
M -t). 
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For isolated neutron stars regions of back-bending exist 
(in Fig. |S] for 2.1M Q ) and are denoted by thick dashed 
lines. The tiny thick-line region of spin-up of the star by 
angular momentum loss just before the onset of the axi- 
symmetric instability (open circle), is a direct consequence 
of stability condition in GR (for other examples of this 
phenomenon, see Figs. 7-9 in Zdunik et al. 2004). 

As one sees, there is no back-bending phenomenon for 
the accreting NSs, even in the most favorable case Xi = 1, 
where the spinning- up of NS is the most efficient. This is 
due to a basic difference in the spin evolution of isolated 
and accreting NSs. For accreting NSs, increase in J is nec- 
essarily accompanied with an increase in Mb . This clearly 
weakens the impact of rotation on the stellar structure 
because of a simultaneous increase in / and in the gravi- 
tational pull. Therefore, in this case spin-up by accretion 
cannot influence the structure of rotating NS so strongly 
as to enforce the back-bending behavior. 

The total baryon mass accreted during the accretion 
spin-up epoch is obviously limited by the mass of the com- 
panion in the LMXB. This results in an upper bound on 
the final frequency reached at the end of accretion. It may 
happen, for example, that the spin-up tracks in Fig. 3 ter- 
minate well before the instability points. 

4. Conclusion 

We studied the spin-up of accreting NSs with hyperonic 
cores and, in particular, the impact of the hyperon soft- 
ening of the EOS of dense matter on the spin-up history. 
The spin-up sequences ended either at the Keplerian (mass 
shedding) limit or, for sufficiently high initial baryon mass, 
by an instability with respect to axi-symmetric perturba- 
tions leading to collapse into a Kerr black-hole. We consid- 
ered an EOS that exhibits a strong back-bending behavior 
in the spin-down of isolated NSs. However, this EOS does 
not show back-bending in the spin-up by accretion, even 
in the most efficient spinning-up, without any angular mo- 
mentum loss in the accretion process. Therefore, for hy- 
peronic EOSs back-bending in the pulsar timing does not 
imply spin clustering in accreting NSs in low-mass X-ray 
binaries (LMXBs). Generally, our calculations show that 
additional increase of mass that accompanies increase of 
angular momentum in accreting NSs, weakens or even sup- 
presses the back-bending phenomenon in the time evolu- 
tion of NS rotation frequency. It should be stressed that 
this suppression of the BB phenomenon for hyperon EOSs 
does not mean that the effect cannot be theoretically ob- 
tained for very strong (perhaps unrealistically strong) soft- 
ening of the EOS. In such a case, the back-bending present 
in the spin-down timing will be also present (albeit in a 
much weaker form) in the accretion spin-up, provided, of 
course, that not too large fraction of the angular momen- 
tum is lost during the accretion process. This aspect of 
the impact of the high-density softening of the EOS on 
the spin evolution of neutron stars and, in particular, its 
observational signatures in the timing of isolated and ac- 



creting NSs will be the subject of our forthcoming publi- 
cation. 
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